Sampling Eulerian orientations of triangular 

lattice graphs 

Paidf Creed 

School of Informatics 
University of Edinburgh 
Edinburgh EH9 3JZ 
Scotland 



Abstract 

We consider the problem of sampling from the uniform distribution on the set of 
Eulerian orientations of subgraphs of the triangular lattice. Although it is known 
that this can be achieved in polynomial time for any graph, the algorithm studied 
here is more natural in the context of planar Eulerian graphs. We analyse the 
mixing time of a Markov chain on the Eulerian orientations of a planar graph 
which moves between orientations by reversing the edges of directed faces. Using 
path coupling and the comparison method we obtain a polynomial upper bound 
on the mixing time of this chain for any solid subgraph of the triangular lattice. 
By considering the conductance of the chain we show that there exist subgraphs 
with holes for which the chain will always take an exponential amount of time to 
converge. Finally, as an additional justification for studying a Markov chain on the 
set of Eulerian orientations of planar graphs, we show that the problem of counting 
Eulerian orientations remains #P-complete when restricted to planar graphs. 

A preliminary version of this work appeared as an extended abstract in the 2nd 
Algorithms and Complexity in Durham workshop. 

Key words: Randomized algorithms, Markov chain Monte Carlo, Rapid mixing, 
Torpid mixing, Eulerian orientations. 



1 Introduction 

Let G = (V, E) be an Eulerian graph, that is, a graph with all vertices of even 
degree. An Eulerian orientation of G is an orientation of the edges of G such 
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that for every vertex v the number of edges oriented towards v is equal to 
the number oriented away from v: degi n {v) = deg out (v). It is well-known that 
the problem of finding an Eulerian orientation of an Eulerian graph can be 
solved efficiently; in this paper we consider the problem of sampling the set 
of Eulerian orientations of a planar Eulerian graph, that is, the problem of 
generating an Eulerian orientation from a distribution that is close to uniform. 
We focus on the Markov chain Monte Carlo method, a standard approach to 
random sampling of combinatorial structures. The Markov chain we study is 
the most natural chain whose state space is the set of Eulerian orientations of 
a planar graph. To move from one Eulerian orientation to another the chain 
randomly selects a face of the graph. If the edges of this face form a directed 
cycle in the original Eulerian orientation the chain reverses the orientation of 
its edges. We will hereafter refer to this chain as the face-reversal chain. 

Markov chain simulation is generally only useful when we know the chain 
is rapidly mixing, that is, when the number of steps required to get within 
variation distance e of the stationary distribution is bounded from above by a 
polynomial in the size of G and e _1 . 

The problems of sampling and counting are closely related; indeed, almost all 
approximate counting algorithms rely on the existence of an efficient sampling 
algorithm. The #P-completeness of counting Eulerian orientations of a graph 
was established by Mihail & Winkler [18], thus motivating the study of rapidly 
mixing Markov chains for this problem. More recently, Felsner & Zickfeld [12] 
have obtained upper and lower bounds on the number of Eulerian orientations 
of any planar map in the more general context of ^-orientations. 

Mihail & Winkler [18] demonstrated that the problem of sampling from the 
Eulerian orientations of any graph can be reduced to the problem of sampling 
a perfect matching of a specially constructed bipartite graph. Sampling a per- 
fect matching of this bipartite graph can be achieved in polynomial time using 
a Markov chain shown to be rapidly mixing by Jerrum & Sinclair [14] (see also 
[15,4]). Nevertheless, studying the face- reversal chain is still of interest as it is 
the more natural approach in the context of planar graphs. In particular, in 
this paper we show that the face-reversal chain provides a more efficient sam- 
pling algorithm for the case of solid subgraphs of the triangular lattice, which 
correspond to configurations of the 20-vertex ice model studied by statistical 
physicists [2]. 

The problem of counting Eulerian orientations of several planar lattices has 
been studied in statistical physics as it corresponds to evaluating the partition 
function Zjce of certain ice models [3,9]. In particular, Baxter [2] has found 
an asymptotic estimate for the number of Eulerian orientations of a grid-like 
section of the triangular lattice. 
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Goldberg et al. [13] and Luby et al. [16] showed that the face-reversal chain is 
rapidly mixing on the set of Eulerian orientations of rectangular sections of the 
square lattice in the cases of fixed and free boundary conditions, respectively. 
Both these proofs followed a similar pattern. First they extended the chain 
with extra transitions. Then a coupling argument was used to find a bound 
on the mixing time of this extended chain. Finally, the comparison technique 
of Diaconis & Saloff-Coste [7] was applied to infer rapid mixing of the face- 
reversal chain. Fehrenbach & Riischendorf [10] studied the mixing time of the 
face-reversal chain on the set of Eulerian orientations of the triangular lattice. 
Following the approach of [16,13] they defined an extension of the face-reversal 
chain, and attempted to use the path coupling technique of [6] to show that it 
is rapidly mixing. However, their analysis does not seem to be correct so this 
problem is still open. 

In this paper we follow a similar approach to [16,13] to show that the face- 
reversal chain is rapidly mixing on the set of Eulerian orientations of solid 
subgraphs of the triangular lattice. In addition to this positive result, we ex- 
hibit an infinite family of subgraphs of the triangular lattice for which this 
chain is not rapidly mixing. Finally, we include a proof that the problem of 
counting Eulerian orientations remains #P-complete when restricted to pla- 
nar graphs as an additional justification for studying the mixing time of the 
face-reversal Markov chain. 

The remainder of the paper is laid out as follows: 

In §2 we define the notation we use for Markov chains and Eulerian orien- 
tations, along with an alternative way to view Eulerian orientations (due to 
Felsner [11]) which proves useful for our analysis. We also describe the ma- 
chinery that we will use to prove results about the mixing time of the chains 
we study. §3 contains a description of two Markov chains, the face-reversal 
Markov chain that is the chief interest of this paper, and an extension of this 
chain which we call the the tower-moves chain. §4 contains the details of the 
path-coupling argument we use to show that this extended chain mixes rapidly 
for solid subgraphs of the triangular lattice. This is then used to infer rapid 
mixing of the face-reversal chain via the comparison method of Diaconis & 
Saloff-Coste [7] in §5. In §6 we show that there exist subgraphs of the trian- 
gular lattice with holes for which the face-reversal chain will always take an 
exponential amount of time to converge. Finally, §7 contains the proof that 
exact counting of Eulerian orientations remains hard when restricted to planar 
graphs. 



2 Background 

In this section we summarise the techniques and theory used in our analysis. 
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A discrete-time Markov chain with transition probability matrix P defined on 
a finite state space Q is called ergodic if it is both aperiodic and irreducible. 
Any ergodic Markov chain has a unique stationary distribution ir. Moreover, if 
P is symmetric then it is uniform over Q. Given two probability distributions, 
p and /i, on Q the variation distance is defined as 

I \p ~ »\ \tv = sup \p(A) - p(A) | . 

The mixing time is a measure of the number of steps taken by a Markov chain 
to get close to its stationary distribution. This is defined as 

r(e) = supminjt : ||P t (x, •) — tt\\tv < e } • 

A Markov chain is said to be rapidly mixing if r(e) is bounded above by some 
polynomial in the size of the elements of Q and in e -1 . For example, in this 
paper we use the number of faces to measure the size of the problem. 

Coupling is a standard technique for proving upper bounds on the mixing 
time of Markov chains [1,6]. A coupling of a chain M. is a stochastic process 
(X t , Y t )t>o on Q x Q such that each of the marginal distributions of (X t )t>o 
and (Y t )t>o are a faithful copy of the original Markov chain. To bound the 
mixing time via coupling we use the coupling inequality [1], which states that 
the variation distance between n and the distribution at time t is bounded 
above by the probability of any coupling coalescing by time t, ie. 

||P* _ n \\ TV < sup p r [X t ^Y t \. 

X ,Yo 

Therefore, in order to obtain a polynomial bound on the mixing time it suffices 
to construct a coupling which will have coalesced (with high probability) after 
a polynomial number of steps. This can be proven by showing that the coupling 
causes all pairs of states to "move together" under some measure of distance. 

We will use a simplified variation of coupling, due to [6], known as path 
coupling. This involves defining a coupling (X t ,Y t ) by considering a path 
X t = Z , Zi, . . . , Z r = Y t between X t and Y t where each pair (Zi,Z i+ i) is 
adjacent in the Markov chain and the path is a shortest path between X and 
Y. This allows us to restrict our attention to the pairs of states which are 
adjacent in the chain, as shown by the following theorem: 

Theorem 1 (Bubley & Dyer [6]) Let M. be an ergodic Markov chain with 
state space Q and let 5 be an integer valued metric defined onQxQ which takes 
values in {0, ... , D}. Let S be a subset ofQxQ such that for all (X, Y) e VtxVL 
there exists a path 

X = Zq, Zi, . . . , Z r — Y 
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between X and Y such that (Zi, Z i+ \) e S for < i < r and 

Y t 8{Z i ,Z i + 1 )=5{X,Y). 

j=0 

Now suppose (X t ,Y t ) is a coupling of M. defined on S. If there exists (5 < 1 
such that for all (X, Y) G S 

E[S(X t+1 ,Y t+1 )\(X t ,Y t ) = (X,Y)} < pS(X t ,Y t ) 

then this coupling can be extended to a coupling (X t ,Y t ) defined on the whole 
of Q x Q such that 

E[6(X t+1 ,Y t+1 )] <(35(X t ,Y t ). 
Moreover, if (3 < 1 then r(e) < log ^^ . 

If (3 = 1 in Theorem 1 then in order to use standard path coupling techniques 
it must be shown that the variance of the distance between any two states 
after one step of the coupling can be bounded away from [6]. However, a 
recent result has removed this condition: 

Theorem 2 (Bordewich & Dyer [5]) Suppose we have a path coupling (X t , Y t ) 
for an ergodic Markov chain M. with distance metric 5 : Q x Q — > N, and 
suppose S C f2 x Q is the set of pairs of states at distance 1. Let t(e) = 
Ip^eD 2 ] ("log e -1 ] , where p is the minimum transition probability between pairs 
of states in S, D is the maximum distance between any pair of states, and let 
t* = Bin(t(e), (1 +p)~ 1 ), where Bin denotes a binomially distributed random 
variable. If f3 < 1 for the coupling, then at the random time t* the Markov 
chain is within e of the stationary distribution, in total variation distance. 

In §4 we will apply Theorem 2 to show that an extension of the face-reversal 
chain is rapidly mixing on any solid subgraph of the triangular lattice. 

To prove a lower bound on the mixing time of a Markov chain we rely on the 
notion of conductance. For any non-empty set S C Q the conductance $(S) 
of S is defined to be 



It is well-known that conductance is inversely related to the mixing time (see 
eg. [20,17]); 

rie' 1 ) > . 

- 2min < 7r (5)<i/2$(5') 

Thus, if n is some measure of the size of the elements of Q then by finding an 

upper bound on the conductance of a Markov chain that is exponentially small 
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in n we can conclude that the chain is not rapidly mixing. Conventionally, a 
chain is said to be torpidly mixing when it is known that it is not rapidly 
mixing. 



Defining the boundary of a set S as 



dS = {x G S : P(x, y) > for some y G Q\S} 



we get 



$(S) < ir(dS)/ir(S) . 



Hence, to show that a chain is torpidly mixing it suffices to find a set for which 
this last expression is bounded above by some exponentially small function. 
This is encapsulated in the following theorem, taken from [17]: 

Theorem 3 If, for some S C Q satisfying < \S\ < \tt\/2, the ratio ix(dS) /ix(S) 
is exponentially small in the size of the elements of Q, then the Markov chain 
is torpidly mixing. 

We will use the above theorem in §6 to show that there exists subgraphs of 
the triangular lattice for which the face-reversal chain is torpidly mixing. 

Let G = (V, E) be a planar graph and F{G) be the set of bounded faces in 
some planar embedding of G. We will use / to denote the number of elements 
in J-(G) and EO(G) to denote the set of Eulerian orientations of G. For any 
simple cycle C C E, Int(C) is defined to be the set of faces on the interior of 
C. A face is said to be directed in an orientation of G if its boundary edges 
form a directed cycle. Felsner [11] showed that it is possible to convert any 
Eulerian orientation of G into another by performing a sequence of reversals of 
the edges of directed faces. Furthermore, a partial order was defined on EO(G) 
by letting X -< Y if X can be obtained from Y by performing a sequence of 
reversals of clockwise directed faces. This order has a unique maximum (resp. 
minimum) element: the unique Eulerian orientation with no clockwise (resp. 
counter-clockwise) cycles. Felsner proved this order forms a finite distributive 
lattice by giving a bijection between the Eulerian orientations of a planar 
graph and a set of functions of the form EO(G) — > N called a-potentials. To 
define these functions Felsner used a partial order on J-{G): k rj if k and 
rj share an edge and that edge is clockwise on k in the minimum orientation. 
The a-potentials are then defined as the set of functions p x : F(G) — > N such 
that: 



The bijection is given as follows: let p x ( K ) equal the number of times k is 



k and 7] share an edge =3- \ p x ( K ) ~ Px(.v)\ ^ 1 
k is on the boundary =^ Px( K ) — 1 

K^jrri=> P x (k) < p x {r]) 



(1) 
(2) 
(3) 



6 



reversed on any shortest path from X min to X . 

In this paper we use Felsner's results to show that the face-reversal Markov 
chain is irreducible and to define a natural distance metric on the set of Eu- 
lerian orientations of a planar graph. For the rest of this paper we will refer 
to this lattice as the Felsner lattice of a planar graph G, and denote it by 
Fels(G). 




Fig. 1. A section of the triangular lattice 

The (infinite) triangular lattice is a structure studied by members of the statis- 
tical physics community, eg. [2,3,9]. This is the infinite graph with the following 
description: 

V = {vij:i,jeZ} 

E = {{vij, v i+i,j}, Hi> v i,j+i}, fad, Ui+ij-i} :i,jeZ}. 

This graph can be embedded in the plane as shown in Figure 1. In this paper 
we focus on two finite subgraphs of this lattice: 

1. A solid subgraph is a graph which can be defined by specifying a cycle 
in the triangular lattice as the boundary and taking everything on its 
interior. 

2. A subgraph with holes is defined similarly to a solid subgraph except we 
do not take everything on the interior of the cycle. 

§4 and §5 are devoted to showing that the face-reversal chain mixes rapidly 
on the set of Eulerian orientations of any Eulerian solid subgraph. In §6 we 
show that this is not true in the case of subgraphs with holes. We prove this 
negative result by exhibiting an infinite family of Eulerian subgraphs with 
holes on which the face-reversal chain is torpidly mixing. 
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3 The Markov chains 

We now define the face-reversal Markov chain on the set of Eulerian orienta- 
tions of any Eulerian planar graph. We use X and X' to denote the states of 
the chain before and after each step. 

One step of the chain M. 

1. With probability 1/2, set X' = X 

2. With the remaining probability 

a. Choose 7 G J-{G) u.a.r. 

b. If 7 is directed then obtain X' from X by reversing the orientation 
of all the edges in 7. 

c. Otherwise, set X' = X 

Irreducibility of this chain follows from the fact that the underlying graph 
of the chain is the cover graph of Fels(G). The holding probability of 1/2 
guarantees aperiodicity, whence the chain is ergodic and converges to a unique 
stationary distribution it. Moreover, the chain is symmetric so tt is the uniform 
distribution on the set of Eulerian orientations of G. 

In order to apply the path coupling theorem we will need extend this chain 
with extra moves in the style of [16,13]. These moves allow us to couple with 
/3 < 1, something which does not seem to be possible for the original chain 
with the metric we use (to be defined in §4). 

Definition 4 Let G be an Eulerian planar graph, X an Eulerian orientation 
of G, and 7 a face of G. We say 7 is almost-directed in X if all but one of 
the edges of '7 have a common direction. We call the edge with the disagreeing 
direction the blocking edge 0/7. For faces k, 77 e F{G) we say there is a tower 
T starting at k and ending at r] in an orientation X if there is a set of adjacent 
faces k — 7 1; . . . , r ) h — rj such that 

• 7j is almost directed in X , and the blocking edge of r ) i is the one shared with 
7»+i for\<i<h-\. 

• r y h is a directed face in X 

We say h is the length of the tower. 

Observe that the definition of a tower implies that C = 1<i</t £(7^) is a 
directed cycle in X . We say that a tower is clockwise (resp. counter-clockwise) 
in an orientation if this cycle is clockwise (resp. counter-clockwise) in the 
orientation. We call r y h and r y 1 the top and bottom of the tower, and refer to 
the right and left sides of the towers in terms of a walk from the bottom to 
the top. It follows that in a clockwise tower the internal edges are all directed 
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from the right to the left, and vice- versa for counter-clockwise towers. 

Let X G EO{G) and let 7 G J-{G) be almost directed in X. If there is a tower 
in X starting at 7 we can find it by walking along the faces of G, starting at 7 
and choosing the face sharing the blocking edge with the current face at each 
step. If at any point we reach a directed face then we have found a tower. If 
we reach an undirected face which is not almost directed, or whose disagreeing 
edge lies on the boundary of the graph, or whose blocking edge is the same 
edge as the blocking edge of the previous face, then there cannot exist a tower 
starting at 7. To see that this process terminates (ie. does not wrap around 
on itself) consider the case where we reach the starting face again. For this to 
occur we must have two cycles C\ and C2, where 

Int(C 2 ) U {7i |1 < % < r} = Int(Ci) and ^ £ Int(C 2 ) , 

with all the edges connecting C\ to C 2 oriented towards C\ or C 2 , something 
which is impossible in an Eulerian orientation. In particular, all towers in an 
Eulerian orientation of a solid subgraph of the triangular lattice are linear 
since any vertex at which a tower bends will be unbalanced. 

We now define the extended tower-moves chain. The definition includes an 
undetermined probability which will be fixed later. 

One step of the Markov chain M. 

1. With probability 1/2, set X' = X 

2. With the remaining probability 

a. Choose 7 G F{G) u.a.r. 

b. If 7 is directed then obtain X' from X by reversing the orientation 
of all the edges in 7. 

c. Otherwise, if there is a length h tower T = (7j)i<i<h with j 1 = 7 
then let C = ®\<i<hFi. With probability pt obtain X' from X by 
reversing all the edges of C . 

d. Otherwise, set X' = X 

This type of chain has been used to extend the face reversal chain in the past, 
see [16,13]. The ergodicity of this chain is inherited from the ergodicity of 
Ai, since every transition in Ai is also a transition in M.. As long as px is 
chosen so that the probability of reversing a tower is independent of whether 
it is a clockwise or a counter-clockwise tower, M. converges to the uniform 
distribution. To see this suppose X can be obtained from Y by reversing a 
clockwise tower T. But then we can obtain Y from X by reversing a counter- 
clockwise tower containing the same faces as T, so P(X, Y) = P(Y,X). 
Hence, the stationary distribution of M, is also the uniform distribution on 
the set of Eulerian orientations of G. 
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4 Rapid Mixing of M on solid subgraphs 

In this section we use the path coupling technique of Bubley & Dyer [6] to 
show that the Markov chain M is rapidly mixing on EO(G) when G is a 
solid subgraph of the triangular lattice and is chosen appropriately. This 
is not the first result regarding the mixing time of this type of chain. M. has 
been shown to be rapidly mixing on the square lattice (using different pr 
values) with different types of boundary conditions in [16,13]. Fehrenbach & 
Riischendorf [10] attempted to give a proof of rapid mixing for a related chain 
(in which only towers of length 2 were used). However, there are errors in the 
analysis of [10], and it seems as though their path coupling analysis cannot be 
fixed for the particular chain used. 

Note that rapid mixing proofs for this chain are dependent on the correct 
choice of the probabilities pt- For example, the proof of [13] sets pt to l/4/i if 
the tower runs along the boundary, and l/2h otherwise, where h is the length 
of the tower T. We now state the main result of this section but defer the 
proof until we have presented some useful lemmas. 

Theorem 5 For any solid subgraph of the triangular lattice the Markov chain 
M. is rapidly mixing with 

T~( e )GO(/ 4 log 6- 1 ), 

when Pt = l/3h for all towers T of length h. 

In our application of path coupling to proving rapid mixing of Ai we will need 
a metric on the set of Eulerian orientations of a planar graph G. We define 
4>{G) to be the covering graph of Fels{G). With this definition, 4>{G) is a 
connected graph with edge set 

{{X, Y} : X can be obtained from Y by reversing a single face} . 

We now define our metric 8 : Q, x Q, — > {0, 1, . . . , D}\ 

S(X , Y) = min{| ^ x ,y I : 1x,y * s a path from X to Y in 4>{G)} . 

In particular 5(X , Y) — 1 if X and Y differ on a single face. Let X min and 
I max denote the unique minimum and maximum elements of Fels(G). The 
fact that Fels(G) is a distributive lattice implies D = 5(X min , X max ), since 
we can always find a path of length D between two states X and Y either by 
going down from X to X min and then up to Y, or by taking a similar path 
through X max . 

Lemma 6 The maximum distance between any pair of Eulerian orientations 
of a solid subgraph of the triangular lattice G with f bounded faces is 0(f?). 
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Fig. 2. Gi, G2, and G3 laid over each other 

PROOF. From the definition of the bijection between Eulerian orientations 
and a-potentials given in $2 we can conclude that the distance between the 
maximum and minimum Eulerian orientations is *)>2~ i eF(G) Pmax.il)- 

(1) and (2) imply that p max (7) is exactly the minimum number of edges a line 
has to cross to going from 7 to the boundary. Let Gk be the smallest graph 
which contains a face 7 with p max (7) = k. We can construct Gk inductively, 
starting with G 1 — K 3 . To extend Gk to Gk+i we add a cycle containing all the 
vertices on the boundary ofGk along with any vertices necessary to ensure that 
Gk is a solid subgraph of the triangular lattice, see Figure 4 for an example of 
how this works. A simple inductive argument shows that the number of faces 
added at each step is 3k. This implies | G Q(k 2 ), so p max (7) G 0(y/J) 

for any face 7. 

In the following let X, Y denote two Eulerian orientations of G which are 
adjacent in <f>(G), let 7 be the face on which they disagree, and let iV(7) 
denote the set of faces which share an edge with 7. Let k be the face chosen 
at some step in the coupling. We then define the coupling as follows: 

• with probability 1/2 — 1/2/, both chains remain unchanged, 

• with probability 1/2/ chain 1 moves using k — 7, and chain 2 remains 
unchanged, 

• with probability 1/2/ chain 2 moves using k — 7, and chain 1 remains 
unchanged, 

• for each k ^ 7, with probability 1/2/ both chains attempt to move using 

K. 

In order to apply the path coupling theorem we need to show that the expected 
distance between the two chains does not increase in a single step of the 
coupling. To do this we need to consider which choice of faces will cause the 
distance to increase, which will leave the distance unchanged, and which will 
cause the distance to decrease. We say the move at k involves 77 G ^(7) if 
k = r] or there is a tower starting at k that contains r\. The distance between 



the coupled chains may increase iff k ^ 7 but the move at k involves some 
77 G iV(7). In the following, we use to denote the maximum total expected 
change to the distance between the states of the coupled chains after a single 
step resulting from moves involving the face 77 in either chain. 

Lemma 7 Suppose rj G N(j) is a directed face in X . Then K[5 V ] = -h. 



PROOF. Let k G F(G) such that selecting k gives a move which involves rj 
in at least one of the coupled chains. We have two cases to consider. 

Case k — T]'. Since rj is a neighbour of '7, and rj is directed in X , it follows 
that the blocking edge of rj in Y is the edge shared with 7. Then T = {77,7} is 
a tower of length 2 in Y with a reversal probability of |. Hence, the coupling 
amounts to reversing rj in X and T in Y with probability | ; and reversing 
rj in X but leaving Y unchanged with the remaining probability. The former 
results in coalescence, whereas the latter yields a pair of orientations which 
are distance 2 apart. Thus, since rj is reversed in X with probability 1, the 
expected value of 5 V is 



1 



(+D(1- |) - K|) 



1 

3? 



Case k, ^ rj'. If the move at k involves rj then there is a tower T starting at k 
in at least one of X and Y which contains rj. Since 7 and rj are both directed 
in X it follows that there must be a tower 7\ starting at k and ending at rj 
in Y that does not contain 7, and a tower T 2 = T\ U {7} starting at k in 
X. Let h be the length of "T 1; so T\ is reversed in X with probability ^ and 



T 2 is reversed in Y with probability 



3(h+l) ■ 



Observe that if we reverse T\ in 



X we obtain an orientation which is distance h + 1 from Y , but if we also 
reverse T 2 in Y then we have the same orientation in both chains. Therefore, 
the expected value of 5 V is 



1 



h i_ L 

\3h 3{h+l)J 3{h+l)_ 



0. 



Thus, the only face whose selection will result in a move which involves rj and 
affects the distance between the coupled chains is rj itself, so E[5 r; ] = 1/3/. 

Lemma 8 If rj E N(j) is not directed in X or Y then M[S V ] is no more than 
_L 

3/" 



PROOF. Observe that, since rj is not directed in X or Y , the blocking edge 
of rj will not be shared with 7 in either copy and that the blocking edge is 
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different in both. Hence, no tower can include both rj and 7. Assuming there 
is a tower containing rj in at least one of the two orientations we have two 
cases to consider. If one of the two orientations has a tower starting at some 
k 7 UiV(7) that contains rj then there will be no tower containing rj in the 
other orientation. This is because the tower construction algorithm described 
in §3 will reach a pair of faces who share the same blocking edge (rj and k). 
Thus, when we are in this situation we can assume that one of the orientations 
will be unchanged after one step of the coupling. In the second case there may 
be a tower starting at rj in either orientation. 




Fig. 3. Example from Lemma 8 

Case 1: Suppose (w.l.o.g.) X is the orientation with a tower containing rj. 
Since no move involving rj is possible in the other orientation we only 
need to bound the total expected change to the distance between X and all 
X' which can be obtained by making a move involving rj in X . 

We begin by showing that any tower containing rj in X must start at 
rj or a neighbour of rj. To see this suppose we have a tower containing p, 
k, and rj where k is a neighbour of rj and p G N(K,)\{rj}. Let u, v, w be 
the vertices of rj, and suppose that the edges of rj are oriented uv, uw, 
and w v. Assume also that rj shares {v, w} with k, and rj shares {u, w} 
with 7. Then p must contain w (recall all towers are linear). To satisfy 
the definition of a tower p must have two edges oriented away from w. 
But this implies that there are 4 edges oriented away from w in X , a 
contradiction (eg. see Figure 8). A similar argument holds for any other 
configuration of the edges of rj. Hence, any tower containing rj must start 
at k, or rj. 

Let h be the length of the tower starting at rj. Then, in the coupling, 
the move in which rj is the chosen face is made with probability • 
and the move in which k is the chosen face (if it exists) is made with 
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probability • 3 r^ + u ■ Since these moves increase the distance by h and 
h + 1, respectively, we have 

_ r _ . h h + 1 1 

Case 2: In the worst case, we could have a tower starting at i] in both orien- 
tations. Let T\ and T 2 denote the towers in X and Y , respectively, and 
let hi and hi be the lengths of each tower. If hi < h 2 then the coupling 
reverses the tower in both chains with probability and only reverses 
the tower in X with probability — ^ . The first of these will give two 
orientations which are distance hi + h 2 + 1 apart, and the second gives 
two orientations which are hi apart. Hence, 



1 

3? 



The analysis is identical if h 2 < hi . 



Combining these two cases, we see that the expected value of 5 V is no more 
than -jj. 



We are now ready to combine the previous lemmas into a proof that the 
Markov chain M. is rapidly mixing. 

PROOF. [Proof of Theorem 5] Consider the coupling defined earlier. With 
probability 1/f one of the two chains reverses 7 causing the two chains to 
coalesce. Combining this fact with the results of Lemmas 7 and 8 we find that 
for all X and Y differing on the orientation of a single face 

E[5(X t+1 ,Y t+1 ) - S(X t ,Y t )\(X t ,Y t ) = (X, Y)} < 3^ - - = . 

It follows that (3 = 1 for the path coupling so we can apply Theorem 2 with 
p = 1/2/ and D E C>(/§) (Lemma 6) to yield r~(e) < Bin(t(e), (1 + p)' 1 ) 
where 

t(c)G O^togc" 1 ). 
But (1 + p)^ 1 — > 1 so, for all practical purposes, we can conclude that r~ e 
O^loge- 1 ). 



5 Comparison of t m and r 



M 



In this section we use the comparison method of Diaconis & Saloff-Coste [7] 
to infer a bound on the mixing time of M. from the bound on the mixing time 
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of M. obtained in the previous section. We will use the formulation of the 
Diaconis & Saloff-Coste result from [19], restated here for convenience. Note 
that we are using E(P) to denote the set of edges corresponding to moves 
between adjacent states in the Markov chain with transition matrix P. 

Theorem 9 ([19, Proposition 4]) Suppose P and P are the transition matrices 
of two reversible Markov chains, M. and M., both with the state space Q and 
stationary distribution it, and let 7r* = min xg n7r(^)- For each pair (u,v) G 
E{P), define a path ^ uv which is a sequence of states u = uo, u±, , . . . , Uk = v 
with (ui,u i+ i) G E(P) for all i. For (x,y) G E(P), let T(x,y) = {(u,v) G 
E(P) : (x,y) G 7u J- Let 

A= max < . r V 1 7 U „ \ir(u)P(u,v) > . 

Suppose that the second largest eigenvalue, \±, of P satisfies \\ > 1/2. Then 
for any < e < 1 

T M {e) G 0{Ar M {e) log 1/tt*) 

To use the comparison method to bound the mixing time of the face-reversal 
chain we need to show that every move of Ai can be simulated by moves of 
the chain M.. Suppose T = {7 1 ,7 2 , • • • , 7^} is a tower in X and that Y is 
the orientation obtained by reversing T. Observe that by the definition of a 
tower 7 h is a directed cycle in X. Then we can perform a move on X in M. 
to obtain a new Eulerian orientation X' in which r y h has been reversed. But 
there is now a tower T' = {71,72, • • • , lh-i} i n X' which can be reversed to 
obtain Y. Repeating this process until we reach Y gives a decomposition of 
the tower move into moves of the chain M.. We begin by bounding the size of 
any tower move: 

Lemma 10 Let G be a solid subgraph of the triangular lattice. Then the max- 
imum length of a tower in an Eulerian orientation of G is 0{yfj). 

PROOF. Let T = (7^ 7 2 , . . . , 7^) denote a tower in any Eulerian orientation 
of G. Assume w.l.o.g that T is a clockwise tower. Recall that all the internal 
edges of a clockwise tower are directed towards the vertices on the left, so any 
Eulerian orientation of G must contain a set of h — 1 edge-disjoint directed 
paths linking the left-side vertices to the right-side vertices. Each of these paths 
must go around the top or the bottom of the tower. But each path that goes 
around the bottom (resp. top) contributes 1 to the distance from f y 1 (resp. •y h ) 
to the boundary. Hence, by (1), (2) and (3), max(p max (7 1 ), p max (7j) > h/2. 
But p max (7) G 0{\ff) for any 7 (see proof of Lemma 6), whence h G 0{yJJ). 
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We are now ready to prove our rapid mixing result for M.. 



Theorem 11 Suppose G is a solid subgraph of the triangular lattice. Then 
the mixing time of the face-reversal Markov chain M. on EO(G) satisfies 

r M (e) eO(f loge' 1 ). 



PROOF. Let P and P denote the transition matrices of M. and Ai. For 

each pair of states (x,y) that differ on the orientation of exactly one face (ie. 
each (x,y) G ker(A4)), we define F(x,y) to be the set of all transitions in M. 
containing the transition t = (x, y) as a sub-move. For each such pair we have 

Ax ' y = it(t)P(t v) £ \luMu)P{u,v) (4) 
= 2f J2 \l uv \P(u,v), (5) 

(u,v)eT(x,y) 

= l + i(|r(x,y)|-l), (6) 

where (5) is due to the fact that all transition probabilities in P are and 
that tt is uniform, and (6) is due to the fact that P(u, v) — if (u, v) = (x, y), 
and P(u, v) = — j- if (u, v) is the reversal of a tower. 

Let k be the face that is reversed in the transition t = (x,y). We need to 
consider the different cases in which t can feature as part of the decomposition 
of a tower move (u, v) G E(P). Observe that there are three different directions 
in which a tower can pass through k and contain t as a sub-move (one for each 
pair of edges of k) . Let 7 and 7' denote the top and bottom of the maximal 
tower (in any Eulerian orientation) which passes through k in the i-th direction 
and whose encoding contains t. We use hi to denote the length of this tower. 
Any other tower which passes through k in this direction and whose encoding 
uses t must be subset of the faces of the maximal tower. Moreover the top ( resp. 
bottom) must lie between 7 (resp. j') and k. Hence, the number of tower moves 
that use a particular transition is 0{tify + 0{h\) + 0{h%) . But hi G 0(y/J) 
whence F(x,y) is in 0(f) for all (x,y) G ker(A^). 

Thus, we have A x>y G 0(f) for all (x,y) G ker(A^). The number of edges 
in G is no more than 3f , so 2 3 ^ provides an upper bound on the number of 
orientations of G, and so also on l/V*. Combining all this with Theorems 5 
and 9 we get 

r M (e)eO(f\oge- 1 ). 
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6 Subgraphs with holes 



Given the small collection of positive results regarding the mixing time of the 
face-reversal chain M. (Theorem 11 and [16,13]) and given that the reduction 
of [18]) allows us to sample from EO(G) in polynomial time for any graph, one 
might hope that M. might be rapidly mixing on the set of Eulerian orientations 
of any planar graph. In fact, this is not true and in this section we exhibit an 
infinite family of subgraphs of the triangular lattice (with holes) for which M. 
is torpidly mixing. Consider the family of graphs H N , of which H 2 is shown 
in Figure 4. Note that the dotted lines in Figure 4 represent the parts of 
the triangular lattice which have been omitted. Formally, is a graph with 
vertex set 

V = {vi : 1 < i < 6(2N + 1)} U {m : 1 < i < 6(2N + 2)} U {wt : 1 < i < 6N} . 
The edges of Hn consists of the disjoint union of three large cycles: 

El = (v 1 , . . . , f 6 (2iV+l), Vl) 
E 2 = («1, . . . , ^6(2^+2), Ui) 
E 3 = « 2 , W U « 3 , V 3 , . . . , U 6( 2N+2), «l) 

eg. see Figure 4. It is the large face in the centre of each of these graphs that 
creates the bottleneck in the Markov chain we will use to show torpid mixing. 
We label this face C and its neighbours r ) i (for 1 < % < QN) . The face that is 
adjacent to both r ) i and 7 i+1 is labelled and the face that is only adjacent 
to 7j is labelled rj^ 

Theorem 12 The face-reversal chain M. is torpidly mixing on EO(H N ) for 
N > 3. 



PROOF. From Theorem 3 we know that Ai is torpidly mixing on a set of 
Eulerian orientations Q if there exists some S C Q, with < \S\ < |^|/2, 
such that \dS\/\S\ is exponentially small in f . 

Let S be the set of all Eulerian orientations X on satisfying px(C) < 1. 
We can define a bijection between S and Q\S by mapping p x to p max — p x 
for each X G S , so \S\ = \£l\/2. 

An Eulerian orientation X is an element of dS if and only if C is a counter- 
clockwise directed cycle in X and px(C) = 1- Eor this to occur we must have 
Pxili) = 1 f or each 7j. Hence, the number of Eulerian orientations satisfying 
this condition is exactly 2 2k since each of the r ) i and i] i can take potential value 
or 1, where k = 2N. 
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Fig. 4. The graph H2 

C is the only directed cycle in X min so \S\ = \S'\ + 1, where S' = {X e : 
Px(C) = 1}. We can partition S' as U/c[fe] Si> where 

Sj = {X en-. Px (C) = l A Pjf ( 7i ) = l 

We can find the size of each of the Si by counting the number of potential 
functions which correspond to members of Si. If X e Si then there are two 
possible values for <p x ( K i) f or eac ^ i with p x {ld = 1 an d Pxili+i) = 1> an ^ 
two possible values for p x (r) i ) f or eac ^ i Pxili) = 1- All of the other K; t 
and r\ i must have potential value 0. Hence, 

\ Sl \ = 2 IW) , 

where c(I) counts the number of circular successions in I. The number of j- 
subsets of {1, ... , k} containing m circular successions is given by the following 
expression: 



c(k,j,m) 



k(j\(k-j-l 
j \mj \j—m—l 



if j — ,j > k , or m < 2j — k 
otherwise 
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Then, 

\S\ = l + Y,^jZ2 m c{k,j,m) (7) 

j=0 m=0 

=£ £ Jfi)G*-t-i)*~ +1+2 " < 8 » 

j=l m=max(0,2j— k) J \ / V / 

>g( 2 /-j 2 ^ (9) 

>(|f^) 2 ^" 3 *f k > 17 ( 10 ) 

> 2^ 2+ T7 )k - 3 (11) 

T/ie /asi (me of this follows from the fact that (pY^J) > 2w when k > 17. 
i/ence ; 

\dS\/\S\ <8-2-w fc 6 0(2^). 



7 Complexity of #PlanarEO 



In this section we demonstrate a polynomial-time reduction from #EO to 
#PlanarEO: 

Name. #EO 

Input. A graph G 

Output. The number of Eulerian orientations of G 

Name. #PlanarEO 

Input. A planar graph G 

Output. The number of Eulerian orientations of G 

This suffices to show that ^PlanarEO is ^P-complete since the #P-completeness 
of #EO is already known [18]. Our reduction uses a recursive gadget and can 
be seen as an application of the so-called Fibonacci method of Vadhan [21]. 

Theorem 13 There exists a polynomial time reduction from #EO to #PlanarEO. 



PROOF. Let G be any non-planar graph for which we have an embedding in 
the plane with I crossings. We start by creating a planar graph G' by replacing 



19 



(a) 



(b) 



Fig. 5. The crossover box 

each crossing {x,y} and {u, v} by a vertex s joined to each of u,v,x,y as in 
Figure 5(a). Using our oracle for #PlanarEO we can count the number of 
Eulerian orientations of G' in polynomial time. Not all of these correspond 
to Eulerian orientations of G, eg. {(u,s),(v,s)} may be present in an Eule- 
rian orientation of G' . We call the configurations of arcs at a crossover that 
correspond to Eulerian orientations of G valid configurations. 

For each k we define Hk recursively as in Figure 5(b) with H given by Fig- 
ure 5(a). Let Gk be the graph obtained by replacing each crossing in G with 
H k , so G = G>. 

Now let Xk (resp. y k ) denote the number of possible configurations of the edges 
of Hk satisfying the Eulerian condition which correspond to valid (resp. in- 
valid ) orientations of G' . These values satisfy 



x k = 4x fe _i + 2y k -i (12) 
y k = Ax k _i + 3y fc _i (13) 

with xo — yo — 1. Now let iVj denote the number of Eulerian orientations 
of G' which have exactly i valid crossover boxes, so iVj = #EO(G). Each 
Eulerian orientation ofG' counted by iVj corresponds to exactly x k y l k l Eulerian 
orientations of Gk, so we can write 

#PlanarEO(G k ) = j^N l x k y l k -\ 

i=0 

Then calculating ^Planar EO(G k ) jy\ corresponds to evaluating the polyno- 
mial p(z) = Y^i=$NiZ % at the point x k /y k . Since x k /y k is a non-repeating se- 
quence (see [21, Lemma 6.2]) it suffices to evaluate ^Planar EO(G k )/y l k for 
k = . . . I + 1 to obtain enough information to recover the values of JVj by 
polynomial interpolation. 

Letting n and m denote the number of vertices and edges in G we have 
\V(G k ) \ = n+(4/c+l)/ and \E(G k )\ = m+(8k+2)l. Since I is certainly bounded 
by a polynomial in n it follows that the whole reduction can be performed in 
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time polynomial in n and m. 



8 Conclusions 



We have shown that the mixing time of the face-reversal chain is 0(f 6 ) on 
the set of Eulerian orientations of any solid subgraph of the triangular lattice. 
This compares favourably with 0(f 7 log 4 (/)), the best known bound on the 
sampling algorithm for general graphs [18,4]. On the other hand we have 
demonstrated that there exist subgraphs of the triangular lattice for which 
the face-reversal chain takes an exponential amount of time to converge. Given 
that the presence of a single very large face is the obstacle to rapid mixing 
in these graphs one might ask whether there exists some function g(f) such 
that the face-reversal chain mixes rapidly whenever the number of edges in 
any face is bounded by g{f). It seems that this question is beyond current 
proof techniques although a proof based on analysing the conductance of cuts 
similar to that used in §6 seems plausible. 

We also showed that the problem of finding an exact value for \EO(G)\ is 
^P-complete when G is a planar graph. This does not preclude a polynomial- 
time algorithm for evaluating \EO{G)\ for more restricted classes of graphs, 
eg. solid subgraphs of the triangular lattice. However, if such an algorithm 
cannot be found then our rapid mixing result could be used to construct 
a fully polynomial randomised approximation scheme, that is, an algorithm 
which can approximate the value of \EO(G) \ to arbitrary precision e and runs 
in time polynomial in / and e -1 . We sketch such an algorithm below: 

Let G be a solid subgraph of the triangular lattice and let 7 be a face 
on the boundary of G. We construct G\ by deleting the edges of 7 from 
G and removing any isolated vertices. It is straightforward to check that 
\EO(Gi)\/\EO(G)\ > 1/4. Hence, we can approximate this value to arbi- 
trary precision with relatively few samples from EO(G). But G\ is also a solid 
subgraph of the triangular lattice so we can repeat this procedure until we 
reach a final G k for which we can easily find \EO(G k )\ ) eg. K 3 . Since the 
number of faces decreases by at least 2 at each step of this process it follows 
that k < f /2. Hence, we can estimate \EO(G) \ to arbitrary precision e in time 
polynomial in / and in e -1 . A detailed description of this general method to 
convert efficient samplers into approximate counting algorithms can be found 
in [8]. 
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